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Abstract 

Using a simple identity between various partial derivatives of the 
energy of the vector model in 0+0 dimensions, we derive explicit re- 
sults for the coefficients of the large N expansion of the model. These 
coefficients are functions in a variable p 2 , which is the expectation 
value of the two point function in the limit N = oo. These functions 
are analytic and have only one (multiple) pole in p 2 . We show to all 
orders that these expressions obey a given general formula. Using this 
formula it is possible to derive the double scaling limit in an alterna- 
tive way. All the results obtained for the double scaling limit agree 
with earlier calculations. 
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1 Introduction 



The O(N) vector model has for some time been an interesting testing ground 
for ideas and techniques that are widely used in the more complicated matrix- 
models. The latter models recently gained importance as theories of two- 
dimensional gravity coupled to matter with central charge c < 1 |]3], |J [J. In 
this article we investigate the large-N expansion of the vectormodel, obtaining 
analytic expressions for the various terms in this expansion. This proces can 
be easily continued to any desired order, leading to exact closed expressions 
in every order. We will show that the results obtained here are in agreement 
with the well-known perturbative expansions of vectormodels, as well as with 
the double- scaling limit, considered earlier in |I[ As an exta bonus the 
method used to obtain these expressions is useful to find non-perturbative 
approximations for the energy of the vector model for any value of N. The 
O(N) vector model is defined by the following partition function: 

E = / d N xe-^ 2 '^ 2 (1) 



e 

A crucial observation is now that for this model all correlation functions 
(or, in this simple case, the moments of the integral) can be obtained by 
successive derivations with respect to the parameter (3. We find for the first 
two correlators: 

--< x 2 > (2) 



dE 

dp 



d 2 E o o , o,o .dE.r, dE 

W =< x y < ( ' >= { 9p ] " ¥ 

We rewrite the last formula in a form that will be of greater use: 

dE = ( dE ]2 _^E 

dg [ d/3 } dp 2 1 ' 

It is indeed this formula that can be used to find a series expansion in the 
coupling constant for any value of N, and will also form the basis for the 
large N expansion of the model. Let us show how this comes about. The 
only required input is the value of the energy for the free model (i.e. without 
the interaction term). Since in this case the integral is simply gaussian, this 
value is easily found to be: E — =j log(/3) + c te . Using this expression as the 
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lowest order approximation to the energy one can calculate with it the RHS 
of (H). We then find in the LHS the derivative of the energy with respect to 
g, given up to zeroth order. Upon integration we find the energy up to first 
order. Explicitly: 



E ^ lom+ m+v g+0i/) ( 4) 



This first order approximation to the energy can again be used in the RHS 
of (|3p, giving us, after integration, the energy to order g 2 . It is clear that 
by just repeating this simple algorithm the energy can be found to arbitary 
large order in the coupling constant for any value of N. The first few orders 
are given by: 



^ log(/3) + N ^ N + 2) g - W + 2){N + 3)j 



N{N + 2)(5iV 2 + 34iV + 60) 3 N(N + 2)(7N 3 + 79N 2 + 3WN + 420) 4 
+ 12^12 9 g^s 9 

(5) 

Since the zero-dimensional integral simply counts the number of diagrams in 
every order of perturbation, we explicitly find the number of diagrams with 
one, two, three, ... index- loops in every order. 



2 Large N expansion 

The simple formula (|^) can however be used to give far more powerfull results 
than just a series expansion in the coupling constant for general N. In fact, 
with a little extra input it is all we need to find explicit expressions for the 
coefficients of the large N expansion of the integral. The large N expansion 
of the vector model consists in taking the number of field components N to 
oo, at the same time rescaling the coupling constant g to g/N. The 1/N 
expansion is then the series expansion around iV = oo of the model thus 
defined. It can equivalently be viewed as an expansion in the number of 
indexloops of a given diagram. For large N the expansion can then be used 
as an approximation to the original integral. It differs significantly from the 
perturbation series in that every order provides us with a nonperturbative 
approximation to the original function. The general formula (|3|) needs a little 
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modification in this case because of the rescaling of g and is easily found to 
be: 

~dg ~ ~N ~d(3 ~ ~d[P { ' 

The following expansion of the energy is proposed: E = NE Q +E% + jjE 2 + 
■■■■ — J^iloN^Ei Putting this into (|6|) and identifying the terms with equal 
powers of we find an infinite series of identities: 

9E _ dE 2 
dg ~ 1 8(3 ' 

dE 1 = dEp dEt d 2 E 
dg 1 d/3 K d/3 } dp 2 

dg { 8(3 )K d/3 ' K d/3 ' d(3 2 1 ' 

and so on, with the p th formula given by: 

9E P _ dEj dEp-j _ d 2 E p _i 
dg ~ ho d P d P d(3 2 U 

We can rewrite this in a more convenient way as: 



9E p _ dE dE p JL^ dEi dE p _i d 2 E p _ 1 

In this form, the second part of the RHS of © consists entirely of lower 
order coefficients in the large N expansion, while only the first term depends 
on E p . It is this property that will allow a recursive determination of the 
coefficients. We know that the energy is a dimensionless quantity, and that 
therefore it can only depend on the dimensionless combination g/(3 2 . We 
thus find the following relation: 

dE P = 2 9dEp (10) 
dp (3 dg 1 1 

Using (H) in © we find: 
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where is defined as being equal to the second term in the RHS of (^) 
and p 2 is the expectation value of jt in the N = oo limit, (i.e. the derivative 
of E with respect to 0). Since p 2 and E can be obtained from the saddle 
point approximation of the integral, this formula will allow us to find all the 
coefficients in the large N expansion using recursion. 



3 Coefficients of the large N expansion 



The saddle point equation of the integral is found by rescaling x to v Nx and 
finding the extrema of the effective action. This gives rise to the following 
equation for p 2 : 

1 - 2(3 p 2 - Agp A = (12) 
for which we choose the regular solution (finite when g — 0): 

P 2 = - ^ + p W 

Taking the derivative of ([12]) with respect to P and g we easily find the 
following identities: 

9p2 p2 pi (14) 



dp (3 + Agp 2 1 - pp 2 
dp 2 -2p 4 -2p 6 



(15) 



dg (3 + Agp 2 1 - Pp 2 

Note that in the last expressions all explicit g-dependence is eliminated. We 
now find an expression for E by integrating the relation p 2 = leading 
to: 

Eo = n^m (16) 

This is all one needs to obtain the coefficients of the large N expansion to 
arbitrarely high order in j^. Using formula (|T5]), we can change integration 
variables from g to p 2 . In doing so we completely eliminate any explicit 
dependence on g. The expressions will be given as functions of p 2 only, (or, 
more correctly, as functions of the dimensionless quantity Pp 2 ). After the 
change of variables formula ([TT|) becomes: 
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24 

The integrals encountered in this problem are relatively simple and can 
all be carried out exactly. Note that this would have been quite more com- 
plicated without the change of variables. Since all lower order coefficients 
are needed to find a coefficient of higher order, we have to determine the 
expressions one by one in rising order. Beginning with E\ we find without 
too much difficulty: 

, P 4 
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= ^ln(2-2/3p 2 ) (18) 

In higher orders the expressions become more elaborate, but no serious com- 
plications arise. The integrals are always of the same type (simple integrals 
of rational functions) and can easily be carried out exactly. For E 2 we find 
in an analogous way: 

1 1 d(3 1 df3 2 
p 4 1 + 2/V - 8(3 2 p A 



2\2 



4 (1-Pp- 
l(2^-l)W + 4) 

E2 = "24 JT^f (19) 

Explicit results for higher order are collected in an appendix. A closer 
look at those results reveals that from E 2 on the coeffcients appear to obey 
some general form: 

^ = ( i%y.) p ^» (2 °» 

Where P p is a polynomial of degree at most 2p — 3. This form holds for all the 
expressions that are given in the appendix. We will now prove by induction 
that it is true for all orders. 
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We already know that the energy is a dimensionless quantity. It thus 
depends only on the dimensonless variable j3p 2 , which we will from now on 
denote as a. The derivative of a function of a with respect to (3 is given by 

— ft ) = df_ffo 
dp f[a) da dp 

= a l~2a df 

(3 1 - a da 1 ' 

If the function / happens to be of the form (PU[), the derivative with 
respect to (3 is given by: 

Of a (2a - If 

dp = P(l-a)^ Qp{a) (22) 
With Q p a polynomial. The second derivative is given by: 

d 2 f a 2 (2a -If 

W = W (1 - a)^ {a) (23) 

Now we assume that all obey the form (|20|), i=l,...,p-l. In that case 
F p _i is given by: 

a 2 (2a -If- 1 

With this form for F p _i, the coefficient E p is given by: 

1 r a (2a' — ~n p_1 
E ? = - 2 h L a0 U ^(»')^ (25) 



It is clear that also this integral is of the form proposed in (p0|) . The integral 
can be performed by expanding the numerator in powers of 1 — a'. The 
integral is then reduced to a sum of simple integrals of inverse powers of (1 — 
a'), the highest power being 3p — 2. After integration the result is therefore 
still a sum of inverse powers of 1 — a, this time with highest power 3(p — 1). 
The integral is is therefore given by a polynomial divided by (1 — a) 3 ^ 1 ^. 
Furthermore, since the integral vanishes for a — ^, the RHS must at least 
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be proportional to (2a — 1). But, since the integrand is proportional to 
(2a' — the integral itself is proportional to (2a — We have thus 

proved by induction that the form proposed for the p th coefficient in the 
large N expansion of the energy of the vector model is indeed correct. This 
general form can be used to obtain the double scaling limit of the vector 
model. Two remarks are at order here: in the proof we assume that Ex 
has the form proposed in (|20|) . This is clearly not the case. However, the 
derivative of E\ has the correct form, and therefore the proof remains valid. 
Secondly, one sees that the general form b striking ressemblance to 

the form conjectured for the coefficients of the large N expansion of matrix 
models. In the following ansatz was proposed for these coefficients: 

with Ph a polynomial and a 2 the solution of a quadratic equation similar to 
(HI. 



4 Double scaling limit 

From ( pOl) we find that all the coefficients exhibit a pole at a = 1. This 
pole corresponds to the branch point singularity of p 2 , namely: g c = —@r-. 
If we define e by e(3 2 = g — g c , a can be developed around its critical value 
and we find: a = 1 — 2^fk + 0(e). Substituting this in the general form we 
have found for the coefficients of the large iV expansion we find that near the 
critical point they behave as: 

E p ~ (ly-'P^e-l^ (27) 

This is of course not true for E and E\. It is however straightforward to 
obtain their expansion around the critical point with the expressions we have 
found: 

E (g)-E (0)= 1 -- l ^ + e- 8 -el + O(e 2 ) 

El = l -\og(A^-e) + 0(e) (28) 
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This is exactly the same behaviour found in |I| (because of a different defini- 
tion of e, one has to replace e in the above formulae with 4e). The behaviour 
of E and E% was dubbed "non-universal" and therefore these terms were 
omitted in the double scaling limit. We will follow this procedure and de- 
fine the energy as E = Y^=2 aFFt-^p- Near the critical point the energy is 
approximated by: 



OO 1 
3,. 1 _ ,^ 



E^J2 N-&'- 1 )(e-5)P- 1 P p (l)(-) 1 '- 1 (29) 

p=2 8 

One remarks that in every term we now have a competition between a sup- 
pression factor proportional with an inverse power of iV and a divergence 
proportional to a negative power of e. The double scaling limit consists in 

3 

taking iV to oo and e to 0, while keeping the combination Ne? finite (this is 
the same combination as found in [l]). In this way we sum the contributions 
from all loops. If we define the double scaling variable as: z = (iVea) -1 the 
energy is given by the series expansion: E — J2 ^ P_1 (|) p ~ 1 -Pp(l)- We will now 
prove that the complete series expansion is indeed the one proposed in [|l|]. 
We find that in the double scaling limit (a ~ 1 — 2^fe) the following relation 
holds: 

df(a) -a df(a) 1 df(a) 
d(3 2^(3 da 2/3 de 1 ' 

Therefore the first and second derivatives with respect to (3 are given by: 

gg, ^ -3(p-l)a §p+ i 
dp ~ 4/3 p ~ l 

g^ 3(p-l)(3p-l)a» 

d/3 2 ~ 16/3 2 p ~ l 1 ' 

(and similarly for Ex, where E p is assumed to have the form b p _x£~^ p ^) 
Another usefull relation is given by: 



which we find by changing variables in fl25|) from a to e and taking the deriva- 
tive of both sides with respect to e. Using these relations in the definition 
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of Fp-x and identifying equal powers of e we readily find a recursion relation 
for the coefficients b p : 

3 3 

- k>\ = tE -?)(?-!)- tp(p - ( 33 ) 

Now, from the definition of the energy in the double scaling limit we find: 

oo oo 

£ = EVi^~ 1 = EV^ (34) 

p=2 p=l 

Therefore the recursion relation (|33|) can be cast into a differential equation 
for the energy. This can be done by multiplying the relation with z p ~ l and 
summing p from 2 to oo. 

oo o oo p— 1 o oo 

- 4 £ p = - E E VA-*(p -?)(*- 1)^ 1 - t E p(p - iJV-i^ 1 

p=2 4 p=3 g=2 4 p=2 



The LHS is immediately found to be equal to: 



(35) 



% + 4 ^i (36) 



While the first term in the RHS is equal to: 



!<fV-.(«-l)««- I ) !, = !<*^)" (37) 



And the second term finally: 



3 z d 2 (zE) 



4 dz 2 

dz 



(38) 



Defining ^ as ^ this leads to the following differential equation: 



At 

3z 2 -^- - 3z 2 £ 2 + 6z£ - 16f = 5/12 (39) 



One can also derive an equation for the partition function ( = e 



-E. 



Z 



d 2 ( _ d( lQd( 5 



a- +22 9l-YS + 36 C = ° (40) 

This is exactly the same equation as found in |IJ , apart from a rescaling from 
z to z/8. We have therefore proved that the double scaling limit defined in 
this way confirms the earlier derivations. 
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5 Conclusions 



In this paper we have explicitly obtained the large N expansion of the vector 
model. The coefficients are found to be analytic functions of a variable p 2 , 
which is the expectation value of the two point function in the limit N = oo. 
Remarkably enough all these functions (apart from the first two) obey some 
general form, which clearly illustrates the divergence of the coefficients. It 
is seen that all coefficients diverge at the same critical value of the coupling, 
and in such a way that a double scaling limit may be taken. This alternative 
derivation is in a sense closer to the philosophy of the double scaling limit 
because it clearly illustrates the competition between the large N suppression 
and the divergence at critical coupling. The results derived in this paper are 
however far more general than their application to the double scaling limit. 
They also allow us to obtain information of the model away from the critical 
region and to find approximations for the energy for finite values of N. 



A Appendix 

We collect some important results in this appendix. The first six coefficients 
in the large N expansion are given by: 

/V - io g (/V) 

E - 

E 1 = Un(2-2(3p 2 ) 

1 (2/V - 1) W + 4) 
2 24 (1 - /3p 2 ) 3 

5 (2/V - 1) W 
3 16 (l-/9p 2 ) 6 
(2/V - l) 4 (8/3 5 p 10 - 56/3 4 p 8 + 164/3 3 p 6 - 4696/3 2 p 4 - 1841/V + 896) 



E, = 



5760(1 - (3p 2 ) 9 
(2/V - l) 5 /V(706/3 2 p 4 + 823/V - 399) 



256(1 - (3p 2 ) 



12 
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The first terms of the series expansions of these expressions are given by: 

, 2 80 „ nn . 1072 _ 
E = 4/ + -/ - 224/ + —/ + ... 

3 5 

2 704 „ „ , 197632 . 
E l = 22-20/ + — /- 2976g 4 + —^—g 5 + .. . 

E 2 = -24/ + —-/- 14976/ + 296960/ + ... 
3 

E 3 = 640/ - 33280/ + 1126400/ + ... 

El = -26880 9 ' + 106 ° 4544 / + ... 

5 

E 5 = 1548288/ + ... 

Notice that the sum J2i Ej gives the correct expansion of the " N — 1 vector- 
model" , which is just a one-dimensional integral: 

9 , 25671168 , 

E = 3g- 48/ + 1584/ - 78336/ + / + ... 

5 

(where we have put f3 — \) 



References 



[i] 

[2] 
[3] 
[4] 
[5] 
[6] 
[7] 



P.Di Vecchia, M.Kato, N.Ohta,Nucl Phys B357 (1991) 495 

J. Zinn- Justin, Phys Lett B257 (1991) 335 

M.Douglas, S.Shenker, Nucl Phys B335 (1990) 635 

E.Brezin,V.Kazakov, Phys Lett B236 (1990) 144 

D.Gross, A.Migdal, Nucl Phys B340 (1990) 333 

P.Di Vecchia, M.Moshe, Phys Lett B300 (1993) 49 

D.Bessis, C.Itzykson, J.B.Zuber, Adv in Appl Math 1 (1980) 109 



12 



